Fermions with attractive interactions on optical lattices and implications for correlated systems 
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In this paper we address the behavior of the superfluid transition temperature T c in the attractive Hubbard 
model. The attractive Hubbard model can be regarded as the generalization of BCS to Bose Einstein conden- 
sation (BEC) crossover to a lattice and as such, may have implications for future optical lattice studies. Never- 
theless, the BEC limit of the Hubbard model is very different from that of the more well studied Fermi gases, 
owing to the strong inter-site repulsion between pairs. Here we address systematically the effects of pairing 
fluctuations for all filling fractions over the entire range of attractive interaction strength. A central conclusion 
of our work is that in a lattice, around half filling, the smooth evolution from the BCS to the BEC limits is in- 
terrupted: T c vanishes when the system approaches the bosonic regime with increasing interaction strength. We 
suggest that this interruption of crossover may signal a quantum critical transition to another form of superfluid 
not continuously connected to a BCS-like phase. A simple variational ansatz for an alternate ground state in 
this more strongly coupled superfluid is presented. A generalization of the (s-wave) Hubbard model to d-wave 
pairing allows us to address issues of relevance to high T c superconductivity. The phase diagram (representing 
the pairing or pseudogap onset temperature T* and T c ) shows that here too, one observes a vanishing of T c 
when T* becomes sufficiently large. Given this predicted breakdown of the crossover, and given the striking 
similarity to features of the cuprate phase diagram, we suggest that future experiments on ultracold fermions 
in optical lattices should not be exclusively limited to the repulsive Hubbard model, but should address the 
attractive model in order to elucidate features of high temperature superconductivity. 

PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 



I. INTRODUCTION 

Recent experiments on ultracold Fermi gases of 40 K and 
6 Li H in the presence of an optical lattice are creating con- 
siderable interest in the community. These experiments can be 
viewed as quantum simulations of important condensed 
matter problems such as the repulsive and attractive fermion 
Hubbard models. In these experiments one can vary (albeit, 
not independently) both the on-site interaction strength be- 
tween fermions and the lattice depth and in this way change 
the interaction U and hopping t. In the case of attractive in- 
teraction U, one can view the Hubbard model as generaliz- 
ing the important problem of BCS- Bose-Einstein condensa- 
tion (BEC) crossover 0, 01 to the case of a lattice. This 
crossover study represents an exciting research area because it 
provides a means of extending what is arguably the paradigm 
of all theories in condensed matter physics (BCS theory) to a 
far more general situation. 

Many have argued that 18J the key to high temperature su- 
perconductivity comes from the repulsive Hubbard model and 
there is a growing impetus to use cold atoms to address the 
issue of whether this repulsive model can ultimately lead to 
d wave superconductivity Jgl [Toll . Nevertheless, our paper 
is based on the premise that we have potentially as much to 
learn about high T c superconductors from the attractive as 
from the repulsive Hubbard case. This attractive model rep- 
resents a generalization of BCS-BEC crossover which is of 
current interest in atomic Fermi systems QjJ, Il2ll to now in- 
clude the presence of an optical lattice. We will, in addition 
to the s-wave case, consider a natural generalization to the d- 



wave pairing symmetry. The former, however, has been more 
extensively addressed in the literature and is of considerable 
value since the experimental counterpart is more straightfor- 
ward. The latter is more relevant to high T c superconductors. 

We note that the BCS-BEC crossover scenario is argued 
lfl3ll to be relevant to high temperature superconductors be- 
cause of their anomalously short coherence length. Additional 
support for the potential relevance of the attractive Hubbard 
model (AHM) to the copper oxides comes from the fact that 
this crossover scenario naturally leads to a "pseudogap" phase 
which is quite prominent in the cuprates, although there is still 
no consensus on its origin. J^, [IJ, [3 The pseudogap state 
corresponds to a regime (above T c ) where there is a gap for ex- 
citing fermionic quasi-particles. Within the crossover scenario 
a natural interpretation of this gap associates it with the exis- 
tence of metastable (sometimes called pre-formed) pairs aris- 
ing from the stronger-than-BCS attractive interaction. These 
noncondensed pairs require an added energy in order to break 
them apart and create fermionic excitations. 

A central goal of this paper is to present a phase diagram 
for the various regions of d- and s-wave superfluid stability 
as a function of interaction strength and band filling. In this 
context, we also present a phase diagram indicating how the 
characteristic temperatures, i.e., the pseudogap onset temper- 
ature T* and the condensation temperature T c , vary with the 
strength of the attractive interaction in the <i-wave case. When 
the attraction gets progressively stronger (as measured by the 
size of T*), T c begins to decrease and ultimately vanishes. 
This behavior is strikingly similar to that found in the cuprate 
phase diagram ([8], especially Fig. 7) where, via larger T*, 
one can associate stronger attraction with reduced doping con- 
centration. These observations suggest that future cold gas 
experiments should not focus solely on the repulsive Hubbard 
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model, but also give extensive attention to the attractive case. 

In so far as they relate to the high temperature supercon- 
ductors, the attractive and repulsive Hubbard models should 
be directly compared. In particular, these two models may be 
essentially equivalent if in the attractive case, the interaction 
is presumed to be d-wave. The repulsive Hubbard model is 
thought to give rise to a d-wave pairing iHHOH] although this 
has yet to be conclusively demonstrated. Given that the pair- 
ing symmetry in the cuprates is known to be d-wave and that 
the pair size is anomalously small, it is certainly of equal inter- 
est to study a d-wave generalization of the attractive Hubbard 
model, as we do here. Within this latter context no assumption 
is made about the microscopic origin of the attractive interac- 
tion. Quite possibly, however, this inter-site attraction comes 
the Coulomb interaction which is repulsive at short distances, 
but can be even more complex than contemplated il5l[l6ll by 
the simple repulsive Hubbard model. 

While there have been a number of numerical studies of 
the s-wave pairing model 01 HHil, alon g wi th approaches 
based on dynamical mean field theory ll20l l2l[ I22I l23ll . ours 
is the only systematic study over all filling fractions and over 
the entire range of attractive interactions. We investigate the 
stability of the BCS-Leggett like ground state (applied to the 
lattice). This state is a natural one to consider for the purposes 
of ultimately shedding light on the cup rates, since there is a 
strong belief in the community H 1 1 lL [Till that when super- 
conductivity is present, it is in many aspects not so different 
from the d-wave BCS state. In this way we argue that, while 
more sophisticated ground state wave functions might be of 
great interest for reaching a general understanding, they may 
not be of specific relevance to the cuprates. 

We find that, even though the strict T = mean-field equa- 
tions suggest that superfluidity can occur everywhere, fluctua- 
tion effects lead to a vanishing T c over an extended range suf- 
ficiently near half filling and for moderately strong attraction. 
An important clue underlying this breakdown of BCS-BEC 
crossover is that it occurs very close to the point where the 
fermionic chemical potential fi changes sign. We naturally 
interpret the /1 = point as corresponding to the transition 
into a regime where the system is effectively bosonic. Here 
the "bosons" are associated with condensed or noncondensed 
(finite momentum) Cooper pairs. 

To gain further insight, we are motivated, then, to con- 
sider an effective bosonic model representing the AHM. In- 
deed, when one expands this Hamiltonian in terms of t 2 /\U\ 
the model which emerges contains bosonic hopping as well 
as an inter-site repulsion of precisely the same magnitude. 
At low density the repulsion is not important, and the sys- 
tem is expected to support BCS-like superfluidity with a T c 
of order t 2 /\U\, as originally conjectured in the important 
paper by Nozieres and Schmitt-Rink (NSR) ||6[]. However, 
closer to half filling in the t 2 / [/-expanded bosonic Hamilto- 
nian the inter-boson interactions are strong. We emphasize 
that in this simple BCS-Leggett ground state the interaction 
between bosons is implicitly taken to be weak. More sophis- 
ticated wave functions are required to capture the effects of 
stronger inter-boson interactions 02411 which must go even be- 
yond Bogoliubov level theory to be consistent with the physics 



of the t 2 /[/-expanded Hubbard model. This, then, provides 
an explanation for the failure of this simple BCS trial wave 
function to support a superfluid phase in the BEC regime on 
a lattice. Finally, these studies are extended to the <i-wave 
case as well, where the breakdown of the BCS-Leggett like 
superfluid occurs well before the chemical potential becomes 
negative. 

There has been a substantial literature on BCS-BEC 
crossover on the subject of Fermi gases OEIH with tun- 
able attractive interactions. There is also considerable inter- 
est in whether an analo gou s crossover can be observed for 
fermions on a lattice |f25t|26ll. To effect this crossover we note 
that there are two ways of increasing the dimensionless attrac- 
tion U/t, either by changing the lattice depth or changing the 
on-site attraction U by exploiting a Feshbach resonance. For 
the present purposes of addressing BCS-BEC crossover in the 
context of the AHM |0,|27t] (which we have argued above may 
contain information about high temperature superconductors) 
we should not exploit these Feshbach effects. Rather, the hop- 
ping t and the on-site attraction U are varied by changing pa- 
rameters of the optical lattice. Thus, our conclusions about the 
interruption of BCS-BEC crossover near the onset of BEC are 
entirely associated with the Hubbard model and do not pertain 
if one uses a Feshbach resonance to create bosonic degrees of 
freedom, prior to applying an optical lattice. In this alterna- 
tive circumstance the more appropriate many body model is 
the Bose Hubbard model which we do not consider here. 

There is an extensive literature which addresses the 
AHM, (lHH and references therein), principally in the 
two-dimensional case because of associations with high- 
temperature superconductivity. In a related fashion, the pos- 
sibility of observing (bosonic-like) Mott insulating states at 
full filling has been the motivation for studies from a num- 
ber of groups ll29l [3^, I3TI [32ll . particularly from the atomic 
physics community. There has also been a focus on charge 
density wave states which are energetically degenerate with 
the superfluid phase at precisely half filling and may compete 
or co-exist with superfluidity l28l [33ll slightly away. Indeed, 
a supersolid phase contemplated in this literature is viewed 
as a mixture of charge density wave and superfluidity. Fi- 
nally, we note that there is also a methodology for address- 
ing the AHM via a mapping to the repulsive Hubbard case. 
Except at exactly 1/2 filling, the counterpart repulsive model 
must be solved in the presence of a rather complicated con- 
straint, which is difficult to implement. A general theorem 
which states that, in the AHM, the ground state contains no 
magnetic order [34], must be imposed B28I1 in any study based 
on this mapping. 

For the most part these previous studies have been at zero 
temperature. Here we focus on general temperatures, T. By 
choosing to address the superfluid transition temperature, T c , 
we effectively introduce fluctuation contributions in a fash- 
ion that is consistent with the BCS-Leggett ground state and 
capable of addressing finite temperatures. Without these fluc- 
tuations, in the s-wave case, it appears that a BCS-Leggett 
superfluid ground state is stable for all parameters. By con- 
sidering the entire range of (attractive) interaction strengths 
and filling factors, we enter into regimes which have also not 
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been addressed by complementary numerical or mean field 
techniques. 

Our theoretical scheme is based on a particular T-matrix 
approximation for the pairing fluctuations which is compat- 
ible with this simplest ground state. Alternative T-matrix 
schemes have been applied to attractive Hubbard models in a 
number of different variations. The transition temperature T c 
was estimated in [ 17] within an approach in which there is no 
pseudogap phase in the normal state. By contrast, here as in 
earlier work j25l[35l[36ll . we use a scheme in which supercon- 
ductivity emerges in the presence of a pseudogap. Alternative 
approaches based on dynamical mean field theory |20ll2ll , |22ll 
in three dimensions (3D) and on quantum Monte Carlo sim- 
ulations in 3D Hal or in 2D IU9I1 have also addressed the 
size of the transition temperature. For intermediate attraction 
strengths, a maximum in the T c curves is found, which some 
have argued IT261 l27ll corresponds to the regime where BCS- 
BEC crossover occurs. However, we stress that this maximum 
appears deep in the fermionic regime, quite far from where the 
fermionic chemical potential fi becomes negative. 

It should be noted that there are also extensive studies of 
the Bose Hubbard model B37ll38ll39ll which are viewed as rel- 
evant to the strongly attractive regimes of the fermion AHM. 
We emphasize here, however, that there is an important dis- 
tinction between the composite boson BEC limit and the Bose 
Hubbard model, due to the different commutation properties 
of the field operators for fermion pairs versus those for true 
bosons. This difference is less important in the very low fill- 
ing regime, but it cannot in general be ignored. Finally, on the 
subject of BCS-BEC crossover for the d-wave case, there are 
earlier studies ll3~5l l40ll in the literature. In Ref. ll40il a simpli- 
fied model addressed the 2D square lattice at T = 0. 

II. GENERAL FORMALISM: BCS-BEC CROSSOVER ON 
A LATTICE 

The AHM Hamiltonian is given by 

jjAHM _ £ t y cU. + ^5>-l) 2 . (1) 

<i,j>,<r i 

Here denotes nearest neighbors, tij denotes the hopping 
coefficient, a =|, J. denotes spins, U < is the on-site at- 
tractive coupling constant, n,; = + 7^, and m a — (y ia Ci a . 
We consider the case where there is an equal population of 
both fermion spin states. The fermion filling factor is n = 
Yli{ n i)/N, where N is the total number of sites on the lat- 
tice. In the single-band tight-binding approximation only the 
case < n < 1 needs to be considered due to particle-hole 
symmetry. 

In order to include d-wave pairing as well, it is convenient 
to consider a generalization of the AHM in a momentum space 
representation. 

H - ,*N = Y,&<L°k* 

ko- 

+ ^2 ^ k ' C k+q/2T C -k+q/21 C -k'+q/21 C k'+q/2r ( 2 ) 
kk'q 



Here £k = £k — /•*> A* i s tne fermion chemical potential, 
N = c ik c k * s me tota l num ber operator, £k is the en- 
ergy in the tight-binding band. The interaction in momen- 
tum space assumes the separable form Vk. k' = C^aVk'/'k'- 
Here U a — U and ip^ = 1 for s-wave and U a = Ud 
and ipu = cos k x — cos k y for d-wave pairing. For the 
s-wave case, we consider isotropic 3D square lattices with 
£k = 2<(3 — cosk x — cosky — cosfc z ), where t, the hop- 
ping integral, serves as the unit of energy. For d-wave pair- 
ing we consider quasi-2D square lattices with £k = 2in (2 — 
cosfca; — cosfcj,) + 2t±(l — cosk z ), where t\\ and t± are the 
in-plane and out-of -plane hopping integrals, respectively. We 
use tn as the unit of energy, and presume, for definiteness, an 
anisotropic ratio t±/t» = 0.01, which is reasonable for the 
cuprate superconductors l3~5ll . 

When the attraction is weak (\U a \/t — > 0), the superfluid 
phase consists of loosely bound pairs, as in a BCS model. As 
\U a \ is progressively increased the pairs become more tightly 
bound and the system crosses over over to a BEC -based de- 
scription in which there are composite bosons confined to a 
lattice. The simplest possible ground state for describing this 
crossover on a lattice is that associated with BCS-Leggett the- 
ory |a]. We will adopt this ground state here: 

\BCS) = Y[{u k + u k 4 T c_ ki )|Q). (3) 

k 

Here the coefficients u^,v^ = (1 ± ^ k /-B k )/2 and the quasi- 
particle dispersion E k = A 2 ^ . This ground state 

presumes that the condensation is complete and is, more 
generally, appropriate only for weakly interacting bosons or 
Cooper pairs. In the single trap experiments the BEC asymp- 
tote is indeed associated with free bosons and one might ex- 
pect this wave function to be a reasonable starting point. On 
the lattice, in the BEC limit, however, the repulsive interaction 
is equal to the effective kinetic energy of the pairs (as will be 
discussed below). Therefore, in a properly self-consistent the- 
ory, one might expect to see a break down of this ansatz for the 
ground state wave function near half filling where the effects 
of inter-pair repulsion are no longer negligible. 



A. Present T-matrix Scheme 

To solve for T c in a fashion consistent with the known 
ground state constraints 0] of the gap and fermion num- 
ber equations, we follow earlier calculations 
based on a T-matrix theory for pairing fluctuations. The T- 
matrix is given by t(K, K' ,Q) = t(Q)(fkifk/, and contains 
two contributions, t(Q) = t sc (Q) + t pg (Q), describing the 
condensed and noncondensed pairs, respectively. Here we 
take K = (io; n ,k), K' = (iw n ',k'), and Q = (tO m ,q) 
as four-vectors. We use uo n and ft m to represent the Mat- 
subara frequencies for fermions and bosons, respectively, and 
J2k = ^TIuj Ek- ft can be shown that if one takes 
t 8C (Q) = -(aI c /T)5(Q), as T -> the T-matrix theory 
is consistent with the constraints in the BCS-Leggett ground 
state. Here A sc is the order parameter which vanishes at T c 
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but is equal to the total gap at T = 0. The contribution from 
the noncondensed pairs, which contains the sum of particle- 
particle scattering ladder diagrams, is 

'" {Q) " TOO)- <4) 

and the pair susceptibility is x(Q) = J2kG(K)Go(Q — 
Here Gq{K) — (iuj n — fk) -1 is the bare Green's 
function and G(K) is its dressed counterpart. 

To define the appropriate dressed Greens functions, we 
adopt the usual T-matrix expression for the fermion self- 
energy given by E(K) = Gq 1 (K) - G~ X (K) = 
J2q t(Q)Go(Q — K) ( Pk-<i/2- This self-energy can be approx- 
imated as E(if) = -A 2 G (-K)tpl, where A 2 = A 2 
and the pseudogap contribution is 



A 2 

pa 



(5) 



The total fermion 

2J2 k G(K),ot 



number equation is given by 



k 



T k 



2/(S k ) ^ 



6c 

Tk 



(6) 



The BEC condition requires that the pairs have zero chem- 
ical potential at and below T c : i~ x (0) = C/^ 1 + x(0) = 0, 



i.e., 

(0) = </ Q +2^ £g ^k = 0. 



T < T c . (7) 



For T < T c , to satisfy this BCS-like gap equation (which nat- 
urally emerges for A within the present T-matrix scheme), A 
must necessarily contain a contribution from both a non-zero 
superfiuid order parameter, A sc , and a contribution associated 
with noncondensed pairs, A pg . In this way and in the weak 
coupling limit, this approach represents a re-interpretation 
of BCS theory which underlines the strong similarity to the 
treatment of ideal gas BEC. The contribution from noncon- 
densed pairs enters through a gap equation, not a number 
equation, however. The transition temperature T c is the tem- 
perature above which A sc vanishes and is determined by solv- 
ing Eqs. (0, ©, and © self-consistently. 



Below T c the T-matrix may be expanded as t 



aafl — £ 2 q 2 + iTq after analytic continuation. This simpli- 
fies the evaluation of Eq. (0. The effect of the ciifi 2 term 
can be neglected except in a narrow regime near half filling, 
where particle-hole symmetry leads to ao — > 0. The con- 
tribution of pairs is dominated by those near the bottom of 
their energy band so the dispersion is further approximated as 
= £ 2< Z 2 / a o = Q 2 /2A/ P , where M p is the effective mass of 
the pairs on a lattice. Note that T q <C £ 2 q 2 in the long wave- 
length limit so that it is set to in our numerical calculations. 



B. The strong attraction limit 

For the case of s-wave pairing, it is useful at this stage to 
study an approximated Hamiltonian in the limit that \U\/t is 



finite but very large so that the h opping term can be treated as 
a perturbation. Following Refs. IHU 14211 and dropping overall 
constants, the Hamiltonian can be rewritten as 



H, 



eff 



E Jb \ b i + E Jnb - 



rib, 



(8) 



<%,j> 



<%,j> 



Here J = 2t 2 /\U\, h = c^c^, b\ = c^c^, and n bi = b\bi. 
The pair operators b. L and b\ are not strictly boson annihila- 
tion and creation operators, since their commutator £>T] = 
1 — m, and {bi,bi} — 0, where m represents the number 
of fermions at site i. The Pauli principle insures that these 
"bosons" are hard core bosons. This effective Hamiltonian is 
equivalent to an XXZ magnetic model with an effective ex- 
ternal field in the z direction; here the average magnetization 
must have a fixed value equal to (n — l)/2 (see [28]). 

It is important to stress that Eq. ([8]) contains a nearest neigh- 
bor inter-site repulsion which is of the same value J as the bo- 
son hopping term. This inter-site repulsion becomes progres- 
sively more important when the fermion filling is close to one 
half. While the fact that the kinetic energy in Eq. <[8J varies 
ast 2 /U is relatively straightforward to understand, the origin 
of the inter-site repulsion is more subtle. This term arises [6] 
from the energy lowering associated with virtual hopping of 
fermions. Clearly the Pauli principle leads to a constrained 
hopping; if a pair has an occupied nearest neighbor site, then 
hopping will be suppressed, thereby raising the energy. In this 
way we see that there is an effective inter-site repulsion be- 
tween the pairs. 

The Hamiltonian of Eq. ([8]) should be contrasted with the 
boson Hubbard model (BHM) which has been widely studied 
ll37ll in the context of the Mott insulator-superfluid transition. 
The BHM corresponds to true bosons on a lattice, where the 
kinetic energy contribution and the on-site repulsion Ub can 
be varied independently: 



H 



BHM 



(9) 



<i,3> 



III. NUMERICAL RESULTS FOR THE s-WAVE CASE 

The system of equations (|5])-(|7]i can be readily solved to 
yield the transition temperature T c for the Leggett-BCS state 
and with variable t/U and filling factors. We, again, stress that 
there is an important difference [6] between the BEC limit in 
a gas and on a lattice. For the latter the transition temperature 
becomes zero at arbitrarily large attractive interactions. This 
is associated with the fact that the effective mass of the pair is 
infinite; pairs cannot hop without an intermediate unbinding 
which becomes prohibitively costly at large attraction. Very 
early on ||6t] it was anticipated that T c would vary as t 2 /\U\ in 
the deep BEC. We emphasize here that the the validity of the 
scaling, T c oc t /\U\, requires minimally that the fermionic 
chemical potential fi be negative, i.e., the system must be in 
the bosonic regime. Indeed, as in earlier work 113511 we find 
this dependence for T c (but only) in the low density limit. 
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Figure 1 : (Color online) T c curve for s-wave pairing on a 3D lattice 
for n — 0.3 (black solid line). The tail is fitted to the functional form 
t 2 /\U\ (red dash line). Inset: T c curve for n = 0.7. 

Figure[T]presents plots of T c as a function of U/t at n = 0.3 
and n — 0.7 (inset). For the first case, at low filling, T c has 
a maximum in the regime where \i > 0; once fi becomes 
negative we find a long tail with the expected t 2 /\U\ depen- 
dence. This general behavior is consistent with earlier work 
ll25l l35l l36ll . It should be noted that, using rather different 
formalisms from that discussed here, others I20I l23ll have re- 
ported this asymptotic t 2 /U behavior, but have never corre- 
lated it with the bosonic regime, in the sense of requiring a 
negative fj,. The plot in the inset emphasizes this point. At 
high n this figure shows that T c vanishes right before the sys- 
tem crosses over into the bosonic regime. Also labeled on the 
plots for both cases is T" lax , where the transition temperature 
reaches a maximum. As noted above, this maximum occurs 
at intermediate values of \U\/t for all filling. As for the uni- 
tary limit in the case of homogeneous Fermi gases fTllfTill . it 
occurs deep within the fermionic regime, where the chemical 
potential is still positive. 

Figure [2] summarizes our phase diagram in the U/t — n 
plane for the AHM with the shaded regimes indicating where 
the calculated T c vanishes. Also indicated is the location of 
rpmax as s h own m Fig. [T] as a function of n. It is relatively 
constant in n, as observed by other groups lfl7l [l8l f20|] . Be- 
cause of particle-hole symmetry at n = 1, here we only need 
to focus on the n < 1 half of the phase diagram. 

We note that the upper boundary of the shaded region mark- 
ing the breakdown of this superfluid phase is consistently near 
the fi = line (below which Cooper pairs start to behave 
like composite bosons), as might have been expected from 
the previous figure. Interestingly enough, this upper bound- 
ary (say, for n w 0.5) is not so far from the predicted values 
for U/t = 35 at which the superfluid-Mott insulator transi- 
tion takes place for the BHM JH]. Note that for the BHM, U 
must be positive in order to stabilize Mott phases which derive 
from strong on-site inter-boson repulsion. Finally, we note 
that particle-hole symmetry near half filling effectively pins 
the chemical potential near the band center and the system 
has difficulty reaching the bosonic limit where the superfluid- 



Figure 2: (Color online) T — » phase diagram associated with the 
BCS-Leggett state for the AHM. The shaded region shows where 
superfluidity breaks down due to strong pairing fluctuations. The red 
line shows the trace of y™ aa: i n T c plots as shown in Fig. [T] The 
black dashed line shows where fi = 0, corresponding to the onset of 
a t 2 /\U\ dependence in T c , associated with low filling. At n — 1, 
there is an exact particle-hole symmetry. 



ity is suppressed. As a result an increasingly more attractive 
interaction is required to arrive at the shaded region. 

A key observation of our theory is that the character of the 
T c curves changes at the point where fi changes sign. An- 
other key observation is that, at large filling factors, the transi- 
tion temperature vanishes near fi — 0, when the system ap- 
proaches the bosonic regime. This vanishing is associated 
with a localization of pairs, that is, a divergence of the pair 
effective mass. Physically this is not unexpected since we 
have seen the effective BEC Hamiltonian contains an inter- 
site repulsion of the same magnitude as the hopping term. It 
is this repulsion which inhibits pair hopping at large filling 
factors. Only in the low density limit can this inter-site repul- 
sion be neglected, thereby leading to the conventional behav- 
ior, T c oc t 2 /\U\. At larger filling, we find that this simple 
BCS-Leggett ground state will not support superfluidity in the 
|J7| > t limit. 

We end this section by noting that while these calculations 
have been based on a theory which is compatible with the 
BCS-Leggett ground state, we believe our general conclusions 
will also apply if one were to address BCS-BEC crossover 
in a Hubbard model using alternative crossover schemes, for 
example, based on the NSR approach J6[]. Below T c , it has 
been shown [44] that the NSR scheme and related extensions 
treat the pairs in a weakly interacting fashion. This picture 
of weakly interacting "bosons" (at the level of Bogoliubov or 
Popov theory) cannot be appropriate for the Hubbard model 
near half filling and in the limit of strong attraction. Conse- 
quently, a proper self-consistent calculation of the transition 
temperature based on this starting point should show signs 
of the inadequacy of the NSR wave function ansatz, possibly 
through a vanishing of T c such as we have found here. 



Figure 3: (Color online) Transition temperature T c as a function of 
Ud/t\\ for d-wave pairing on a quasi-2D lattice for n — 0.05. At this 
low filling superfluidity can be sustained for some range of attractive 
interactions, even after /i becomes negative at moderately strong cou- 
pling. For strong enough attraction we find a strictly zero T c . Inset: 
T c curve for n = 0.1, for which T c is nonzero only in the fermionic 
regime. 



IV. NUMERICAL RESULTS FOR d-WAVE PAIRING 

In this section we study the transition temperatures for the 
case of d-wave pairing. We can anticipate that as the system 
approaches stronger coupling, pair hopping will be greatly 
suppressed, just as for the s-wave case, and superfluidity will 
shut down. Indeed, because the d-wave pairs are extended 
over two lattice sites, the effects we saw for the s-wave case 
should be even more dramatic. Moreover, because there is no 
standard d-wave analogue of Eq. (8), it is not clear a priori 
whether the fi — point should reflect a qualitative change in 
the physics. 

Figure [3] presents results for the d-wave transition temper- 
ature in a quasi-2D lattice for n = 0.05 and n = 0.1 (in- 
set). A new feature emerges which is not present in the s- 
wave case. In the ground state, there is a threshold in \Ud\/t 
above which weak coupling superfluidity is stable and below 
which it will not survive [40]. Just as for the s-wave case, 
at all fillings there is a maximum in the T c curves which we 
refer to as T™ ax , inside the fermionic regime. Importantly, 
only at extremely low filling (n < 0.1) is T c finite when fi be- 
comes negative. Otherwise T c vanishes for sufficiently strong 
attraction, but still within the fermionic regime, fi > 0. In 
fact, even for those extremely low values of n, where one can 
pass into the bosonic regime, we find that T c reaches zero for 
sufficiently large \Ud\/t. This shut-down of superfluidity is 
associated with a divergence in the pair mass, that is, with 
localization of the pairs, just as for the s-wave case. 

Shown in Fig.|4]is the phase diagram in the Ud~n plane for 
the d-wave case. In the light (gray) shaded region, there is no 
superfluid ground state associated with the BCS-Leggett T = 
equations [5]. This is on the weak coupling side of the phase 
diagram. In the dark (cyan) shaded region, pairing fluctuations 
destroy superfluidity. This corresponds to the strong coupling 



Figure 4: (Color online) T —> phase diagram for fermions on a 
quasi-2D lattice with d-wave pairing symmetry. The lower (cyan) 
shaded regimes shows where superfluidity does not survive as T ~ » 
0. The black dotted line indicates the threshold beyond which super- 
fluidity becomes observable (i.e., in the unshaded regions) and the 
upper (gray) shaded regime corresponds to a normal phase. The red 
solid line shows the trace of T™ ax shown in Fig. [3] The black dash 
line indicates where fi = 0. 



side of the phase diagram, and arises because of inhibition of 
hopping. 

It should be clear that in the d-wave case the regime where 
pairing fluctuations destroy superfluidity is much larger than 
its counterpart in the s-wave case. This reflects the fact that 
d-wave pairing involves nearest-neighbor sites. The finite size 
of the pairs (which cannot be less than the lattice constant) 
distinguishes them from the point-like composite bosons (of 
the s-wave case) and allows them to break and recombine 
with neighboring sites in complicated ways. In the very di- 
lute limit we understand the behavior of Fig. H] as follows. 
a two-body-like bound pair state exists (i.e., when /i < 0), 
the binding energy will be given by E\> w — 2/j, ~ ti 2 /2m£ i 2 , 
where £ is the pair size. Note that in a many-body system, the 
dimensionless quantity E^/Ep must be large for a system to 
be in the BEC limit. The fact that the d-wave pair size £ can- 
not be smaller than the lattice constant sets an upper bound for 
the binding energy. 

Only in the dilute limit where the mean inter-particle dis- 
tance becomes substantially larger than the d-wave pair size 
can one reach sufficiently large values of Eb/Ep to achieve 
a d-wave bosonic superfluid, principally because of the small 
value of n or Ep. In this low density regime, if we continue to 
increase the attraction (at very small fixed n) we again rapidly 
destabilize the superfluid phase. This is associated with the 
fact that increased attraction requires a larger and larger Ej, 
(or equivalently smaller £), which eventually hits the limit set 
by the lattice constant. 

Since the inter-particle distance 1/kp in quasi-2D scales 
as \ j\fE~p ~ nr 1 / 2 , with increasing n (at fixed attractive 
interaction) the relative binding energy Eb/Ep decreases very 
rapidly. As a result, 1/kp soon becomes comparable to the 
pair size £, so that many-body effects strongly suppress the 
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mobility of the pairs. This eventually destroys the superfluid 
state. 

Finally, as the density approaches half filling, additional ef- 
fects arise from particle-hole symmetry. These tend to pin fi 
to its noninteracting value (w 4tii). This explains why the 
chemical potential is high when T c vanishes at densities close 
to half filling. In summary, despite the fact that the mean- 
field ground state equations have a superfluid solution, we 
find that with fluctuation effects included, this d-wave super- 
fluidity (associated with the BCS-Leggett ground state) has 
difficulty becoming established once the attractive interaction 
becomes moderately strong. 

V. LOCAL PAIR STATE (LPS) WAVE FUNCTION AND 
BCS-LPS TRANSITION 

A. The breakdown of the BCS-Leggett state in the \ U\~3> t 
limit 

We have emphasized that the shaded region in Fig.|2]is as- 
sociated with the breakdown of the BCS-Leggett form of su- 
perfluidity. It is natural to ask what is the nature of the phase 
inside the shaded region. We have in the past 112511 character- 
ized this state as "insulating" on the basis of the general "rule" 
that a bosonic system which is non-superfluid is generally lo- 
calized. Moreover, the onset of the shaded region that we find 
marks the onset of an infinite pair mass which is consistent 
with localized bosons. It has been argued 114511 that an insulat- 
ing phase for both fermions and bosons, if it corresponds to a 
finite excitation gap, is only possible when the filling is com- 
mensurate. One way to get around these arguments (which 
allow insulating phases away from commensurability) is to in- 
troduce phase separation which we do not contemplate here. 
It may also be that localized pair states with a "localization 
gap" arising from many body effects, (not disorder), cannot be 
ruled out. Finally, one may also consider other forms of super- 
fluid phases which may be more stable than the BCS-Leggett 
form, although they need not evolve continuously from the 
BCS state. Indeed, based on a study of the true Bose case 
ll39ll there is a suggestion that the ground state might well be 
a superfluid whenever the nearest-neighbor repulsive term is 
smaller than twice the hopping coefficient. 

We observe that when the on-site attraction is strong, 
fermions are expected to form local pairs and one might an- 
ticipate that a better ground state for the superfluid phase is 
one where a given site either has precisely zero (n, = 0) or 
precisely 2 fermions. Singly occupied sites are unfavorable, 
although they do serve as opportunities for virtual hopping. It 
is convenient in what follows to count the number of bosons 
which correspond to 1/2 of the fermion number. Due to the 
particle-hole symmetry about rib = 1/2, we restrict consider- 
ation toO < lib < (1/2) where rib = Y?,i( n bi)/N is the filling 
factor of pairs. 

Following earlier work [I28II . we contemplate a new ground 
state 

|LPS) = Y[(VT^T b + y/nfrt)\0). (10) 



An analogous wave function was discussed for bosonic sys- 
tems 113 811 . We stress that that the correct pair commutation 
relation = 1 — 2n& ( should be used. The ground 

state energy for any < n b < (1/2) is (LPSl^e/zlLPS) = 
Jznb(2nb — l)N. Here z is the number of nearest neigh- 
bors. Following ll38ll the zero-momentum-pair fraction is 
(LPS\bl\LPS) 2 /n b = 1 - n b . Thus, when n b -> nearly 
all pairs are in the zero-momentum condensate. 

In this LPS state, near half filling, roughly half of the pairs 
at T = have finite momentum. Note that this many body 
wave function is associated with a macroscopic occupation of 
the lowest effective single particle energy level, but that the 
effective single particle levels need not be eigenstates of the 
non-interacting system. Thus the ground state in question (for 
the strongly correlated case) is not associated with (k, — k) 
pairing. To see this in more detail, we rewrite the LPS state in 
momentum space 

\LPS) = Y[(VT^Tb + ^4r c u)|0) (ID 

i 

= n(^ /T ^+v^E44i e " l(p+q),Ri )io>- 

i p,q 

Here R, denotes the position vector of the i-th site. From 
this expression one sees that when rib is finite, both finite mo- 
mentum Cooper pairs cj c |C^ k+q | and higher order terms such 

as c ki c k 2 c k c k 4 are important. The presence of these finite 
momentum condensed pairs may relate to the general issue 
of condensate fragmentation [46] which occurs in the pres- 
ence of degenerate ground states. The degeneracy of the strict 
atomic limit |C/|/t — ^ 00 is partially lifted when weak tunnel- 
ing is included, but it appears that a more natural description 
of the superfluid phase should be one which abandons the (k, 
— k) pairing of the BCS-Leggett phase. 

We note that at low filling, one might expect that since the 
nearest-neighbor repulsion is negligible this LPS ground state 
may not be very different from the usual BCS-Leggett state. 
The equivalence is straightforward to establish. For rib *C 1, 

\LPS) = Y[(VT^~b + V^~b4)\0) 

i 

« (l + V^E c kT c -k4)|0> 

k 

« n(^ /T ^+V^cj c T c -ki)l >- d 2 ) 

k 

Here £i b\ = £\ cj^ = £ k 4 T c -ki- when n b -» 0, the 
BCS number equation n = 2rib = 2 ^ k u k and the normal- 
ization condition it k + w k = 1 lead to \/l — rib and 
Uk w y/rib- As expected, the states \LPS) and \BCS) are 
equivalent in the \U\ t provided we consider the n — > 
limit. This is consistent with our earlier numerical results 
which show that there is no breakdown of superfluidity at 
any strong attraction in the BCS-Leggett phase provided the 
filling is low. For small n, the BCS wave function captures 
the main features of the AHM in the two limits \U\/t <C 1 
and \U\/t ^> 1 and it should be appropriate for describ- 
ing the crossover behavior of the AHM for any \U\/t. This 
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Figure 5: Phase diagram for d-wave pairing in a quasi-2D square 
lattice at n = 0.9 (near half filling). Pairs emerge below T* but 
superfluidity only exists in the shaded regime. While the horizon- 
tal axis represents the strength of the attractive interaction, one may 
make contact with the cuprate phase diagram as shown in Fig. 7 of 
Ref. [8|] by noting that the more underdoped the system the larger is 
T*. Thus, one may view the hole doping concentration x = 1 — n as 
a means of parameterizing the pairing attraction. In this way, there 
are some important similarities with the cuprates. 



is no longer the case near half filling in the strongly attrac- 
tive regime. The BCS-Leggett wavefunction does not cap- 
ture the physics associated with strong inter-boson interac- 
tions, thereby leading us to contemplate an alternate super- 
fluid phase, such as that associated with the LPS wavefunc- 
tion. One might speculate that since the BCS-Leggett ground 
state and the LPS ground state describe different types of su- 
perfluids for very different limits, a quantum phase transition 
may occur when \U\/t is tuned between these two limits. We 
speculate that the onset of this quantum phase transition can 
be loosely associated with the boundary curve for the shaded 
region in Fig. [2] 



VI. EXPERIMENTAL IMPLICATIONS 

A. Implications for High Temperature Superconductors 

While there is an extensive literature [8, 9, 10] which ar- 
gues that high T c superconductivity can be "explained" by the 
repulsive Hubbard model, evidence that supports this point of 
view mainly comes from quantum Monte Carlo simulation or 
numerical study using a density matrix renormalization group 
method on a very small (e.g., 8x16) two dimensional lattice. 
Here we want to address what we can learn about the cuprates 
from the AHM. The phase diagram associated with the attrac- 
tive Hubbard case, importantly, extended to d-wave pairing, is 
shown in Fig.[5]fbr the case of n = 0.9, which is appropriate 
for the cuprates and corresponds to hole doping concentration 



x = 1 — n = 0.1. The behavior of the transition temperature 
was discussed in Sec. [TV] Here we add an additional plot of 
T* which corresponds to the pairing onset temperature. This 
is estimated by solving the standard mean field equations at 
A = 0. 

We stress at the outset that the BCS-BEC crossover sce- 
nario neither requires nor contains detailed microscopic infor- 
mation about the pairing mechanism. (More quantitative fits 
of this phase diagram to that of the cuprates have been pre- 
sented in the literature 114711 ). Nevertheless, this phase diagram 
can be generically compared to that of the cuprates (see, e.g., 
Fig. 7 of Ref. |8[]) without any detailed fits when we exploit 
the fact that (in the cuprate data) as x decreases, T* increases. 
The more underdoped the system the larger is T*. Thus, to 
make progress we may view the hole doping concentration x 
as a means of tuning the size of the pairing attraction. When 
T* is matched to the experimental pseudogap onset tempera- 
ture, the d-wave transition temperature displays a maximum 
as in the inset of Fig. [3] just as seen experimentally at the op- 
timal doping. Importantly, one sees that when the attraction 
gets sufficiently strong superfluidity is shut down. Finally, we 
note that at the lower values for T* (that is, in the overdoped 
regime), the superfluid phase in the cuprates is well described 
by d-wave BCS theory and the pseudogap is negligible. Thus 
the end-point for superconductivity at the lower critical dop- 
ing in the cuprates may well be related to the break-down of 
BCS-BEC crossover which we have been discussing in this 
paper. It should be stressed, however, that this disappearance 
of the d-wave superfluid phase occurs well before one reaches 
the BEC limit. Indeed Fig. [4] shows clearly that the fermionic 
chemical potential is clearly positive (and large) when T c van- 
ishes at high fermion density. This is an important issue vis 
a vis the cuprates, since it is clear that these systems have a 
large positive fermionic chemical potential fi and are thus far 
from the bosonic regime. 

A very important issue in high temperature superconduc- 
tors, which has not been resolved experimentally, is a determi- 
nation of the nature of that phase between a vanishing T c and 
a large T* at very low doping concentrations where supercon- 
ductivity disappears. All that is certain is that it corresponds 
to an insulating state J^, [Tl|] in the sense that the resistivity 
appears to increase with decreasing temperature. It should, 
in summary be clear that, although there may be differences 
associated with long range Coulomb effects which may, for 
example lead to phase separation l48ll . studies of attractive in- 
teractions on optical lattices have the potential for elucidating 
important issues in the high T c cuprates. 



B. Experiments on Optical Lattices 

In this section we discuss possible experiments to observe 
BEC-BCS crossover on lattices in the specific context of the 
AHM. It is not at all straightforward to set up this fermionic 
one-band model and to avoid introducing either multi-band ef- 
fects or direct boson-boson hopping. It should be stressed that 
on lattices, BEC-BCS crossover need not rely on a Feshbach 
resonance tuning of the attractive interaction alone. Rather, to 
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simulate the AHM, as shown in Refs. 011! Hi, the on-site 
attraction U and the hopping coefficient t are simultaneously 
tuned and are, thus, not independent. They can generally be 
controlled by varying the lattice potential depth Vb as well as 
the s-wave scattering length a s (of a Fermi gas in the absence 
of lattice potentials). It is conventional to define the recoil en- 
ergy as Eji = h 2 / (2m\ 2 ), where A is the wavelength of the 
laser used to generate the lattice potential and the lattice spac- 
ing is d = A/2. Importantly, when Vb Er and d 3> \a s \, 
it can be shown that the one-band Hubbard model emerges 
with t = (2/y/^)E R £ 3 e- 2 t 2 and U = irE R (a s /d)£ 3 , where 
£ = (V /E R )V* 

In the regime a s > d both multi-band effects and addi- 
tional terms in the Hamiltonian such as density-assisted hop- 
ping terms become important [23J. Here the Hamiltonian can 
no longer be associated with the AHM. BCS-BEC crossover 
can readily occur in this limit without the interruption we have 
found in our Hubbard calculations. Indeed, when a Fermi 
gas is in the strongly attractive regime (the BEC side of a 
Feshbach resonance) in the absence of any lattice potential, 
fermions form pairs with true binding energies (~ ti 2 /2ma 2 ), 
and pairs repel each other with an effective scattering length 
0.612,5 [50]. These two-body effects are robust when an op- 
tical lattice is then "switched on". In this way a BEC limit 
can be readily established. More specifically, one starts with 
a Fermi gas in the deep BEC regime (satisfying a s <C d when 
the lattice potential is turned off), and then increases the lat- 
tice depth. Importantly, the system can be described by a 
Bose Hubbard model with infinite on-site repulsion. In this 
case the two-body binding energy will dominate the on-site 
attractive potential of fermions. There is no effective nearest- 
neighbor repulsion (found in the BEC of fermion pairs) be- 
cause the true binding energy of pairs prevents the presence of 
unpaired fermions and therefore eliminates these virtual pro- 
cesses. It should be noted that the first generation experiments 
from MIT on optical lattices |2[] were performed making use 
of a Feshbach resonance so that before the optical lattice was 
established, the system was near unitarity. Again, these exper- 
iments should not be considered as simulating the one-band 
AHM. Rather they pertain to a model Hamiltonian which is 
clearly different from the fermion Hubbard model. We have 
argued in this paper, as claimed elsewhere 12811 . that this at- 
tractive fermion Hubbard model may have relevance to high 
temperature superconductors. 

From the present perspective, based on the AHM, there are 
a variety of important experiments yet to be done. Among 
these is to see if one can find the predicted breakdown of BCS- 
BEC crossover which may be associated with the transition to 
a new type of superfluid (The BCS-Leggett state to the LPS 
state) or to some other type of non-superfluid order. To do this 
we propose that the Fermi gas is first prepared in the weakly 
attractive regime (on the BCS side of resonance) satisfying 
\a s \ <C d and then the lattice potential is gradually increased 
to suppress hopping. As shown in Fig. 4 of Ref. [23], this is 
equivalent to increasing \U\/t in the AHM when a s < 0. As 
\U\/t increases, a maximum in T c should first be observed. 
Near half filling when \U\Jt is sufficiently large (as predicted 
by Fig. [2j a new type of ordering may appear or alternatively 



(if the associated transition temperatures are sufficiently low) 
the system will be driven normal. It is quite likely, although 
it may be hard to verify, that the loss of superfluid order is to 
be associated with a quantum phase transition near half filling 
and when \U\ » t. 



VII. CONCLUSIONS 

In this paper we have addressed a generalization of BCS- 
BEC crossover in atomic Fermi systems to include the pres- 
ence of an optical lattice. Our specific interest is to address the 
AHM, and in the longer term, its d-wave generalization which 
have been argued to be relevant to high T c superconductors. 
While there have been a number of numerical studies iTTtL fTsIl 
of the attractive Hubbard Hamiltonian, along with approaches 
based on dynamical mean field theory [20, 23], ours is a sys- 
tematic study over all filling fractions and over the entire range 
of attractive interactions. 

We have investigated the stability of the simplest type of 
superfluid phases which represent a natural generalization of 
BCS-BEC crossover which has been widely studied in the 
gas phases. We refer to the ground state wave function as 
the "BCS-Leggett" state and find that, even though the strict 
T = mean-field equations suggest that superfluidity can 
occur everywhere, fluctuation effects lead to a vanishing T c 
over an extended range near half filling and for moderately 
strong attraction. Thus, there is an interruption of BCS- 
BEC crossover which we explain from an analytic viewpoint. 
This interruption occurs very close to the regime where the 
fermionic chemical potential /i = 0. It is a result of the sys- 
tem passing over into the regime where the standard t 2 /U- 
expanded bosonic Hamiltonian (as an approximation to the 
AHM) becomes valid. In contrast with the weakly interacting 
bosons of the BCS-Leggett state (or of the NSR theory), here 
the bosonic degrees of freedom experience a strong inter-site 
repulsion. Importantly superfluidity in these specific (weak 
coupling) forms cannot be supported. 

We posit an alternative superfluid phase which may be more 
appropriate at very strong attraction. These studies have been 
extended to the <i-wave case as well. We stress that the result- 
ing phase diagram has many features in common with that of 
the cuprates. In particular, at sufficiently strong attraction (as 
represented by a large value for the pairing onset temperature, 
called T*), superfluidity is interrupted. This interruption of 
BCS-BEC crossover occurs while the system still has a rather 
large positive fermionic chemical potential, and hence is far 
from the BEC. 

A central goal of this body of work is to make a case for 
future optical lattice experiments to simulate the AHM, first 
for the s-wave and ultimately for the c?-wave case. It is our 
contention that in this way we have as much to learn about 
the cuprate superconductors, as from studies of the repulsive 
Hubbard model which has been conjectured to give rise to d- 
wave attraction. 

This work is supported by NSF PHY-0555325 and NSF- 
MRSEC Grant No. DMR-02 13745. We thank Cheng Chin 
for stimulating discussions and A. Paramekanti for helpful in- 
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